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ZEROES OF PARTIAL SUMS OF THE ZETA-FUNCTION 


D.J. PLATT AND T.S. TRUDGIAN 


Abstract. This article considers the positive integers N for which Cn{^) = 
has zeroes in the half-plane 9?(s) > 1. Building on earlier results, 
we show that there are no zeroes for 1 < A < 18 and for N = 20, 21, 28. For 
all other N there are infinitely many such zeroes. 


1. Introduction 

The Riemann zeta-function is defined as ({s) = for 5R(s) > 1. Through¬ 

out this article we write the complex variable s as s = a + it with a and t real, 
and consider A to be a natural number. Truncation of the zeta-function gives the 

partial sum Ca(s) = 1-1- H-h One may study these partial sums in the 

hope of deducing some information about (’(s)- For a comprehensive treatment of 
these ideas we refer the reader to [5] and [6]. 

Turan [16] showed that the Riemann hypothesis would follow if for all N suffi¬ 
ciently large CAf(s) had no zero in cr > 1. Let be the supremum over all values 
of a for which Cn{s) = 0. Montgomery |9] showed that for all N sufhciently large, 


V'N = 1 + 


(yi-„(i)) 


log log N 
log A 


where the constant I/tt — 1 is best possible. Therefore for A sufficiently large, CAf(s) 
has zeroes in cr > 1. 

Monach |8] made this explicit: for all A > 30 there are zeroes in cr > 1. His 
proof was in two parts: an analytic argument for A > 549, 798 and a computational 
proof for 30 < A < 549,798. The latter proof is contained in [5J Lem. 3.14, pp. 
134-135]. Monach’s work can be combined with the results of Turan and Spira to 
give the following table. 


Table 1. Zeroes of Ca(s) in a > 1 for various values of A. 


Range of A 

Are there zeroes of Cn{s) in lR(s) >1? 

1-5 

No, dH pp. 7-8] 

6-9 

No, [13 P- 550] and [13 Table II, §4] 

19 

Yes, [13 Table III, §4] 

22-27 

Yes, [13 Table III, §4] 

29-50 

Yes, [13 Table III, §4] 

> 51 

Yes, [3 Thm. 3.8] 


We are grateful to Andrew Booker with whom we discussed several aspects of this problem 
and to the referees for their careful reading of our manuscript. The second author was supported 
by Australian Research Council DECRA Grant DE120100173. 
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Indeed, van de Lune and te Riele m actually computed some zeroes of Cat (s) for 
N = 19,22 — 27, 29 — 35, 37 — 41,47. Adapting Bohr’s theorems on values assumed 
by Dirichlet series, Spira [Til Thm. 3] (see also, [TH p. 163]) showed that if Cn{s) 
has one zero in cr > 1 then it has infinitely many such zeroes. 

Therefore, all that remains is to investigate whether, for, 

(1) N e {10,11,12,13,14,15,16,17,18,20,21,28}, 

CAr(s) has zeroes in cr > 1. We find that there are no zeroes for each of these values 
of N. Combining this with Tabled] one proves the following theorem. 

Theorem 1.1. For 1 < TV < 18 and N = 20,21,28 there are no zeroes of(is[{s) in 
the region a > 1; for all other positive N there are infinitely many such zeroes. 

2. Numerical Computation 

2.1. Interval Arithmetic. Almost all real numbers are not exactly representable 
by any finite-precision, floating-point system such as the 64 bit IEEE implementa¬ 
tion available on most modern processors. Thus any computation involving such a 
floating point system will, unless we are very lucky, only produce an approximation 
to the true result. One way of managing this is to use interval arithmetic (see, for 
example, [10] for a good introduction). Instead of storing a floating point number 
that is an approximation to the value we want, we store an interval bracketed by 
two floating point numbers that contains the true value. 

Interval arithmetic has been used to manage the accumulation of round-off and 
truncation errors. In this paper, we exploit the technique to get zero free regions 
rigorously. As an example, consider the function / : R —> K defined by 

f{x) = — Ax -\- Z. 

Suppose we wish to demonstrate that / has no zeroes for x G [4, 5]. Then we can 
compute 

/([4,5]) = [16,25] - [16,20] + 3 = [-1,12]. 

Since this is inconclusive, we try again, but this time with the interval split in two. 
We have 

/([4,4.5]) = [16,20.25] - [16,18] + 3 = [1, 7.25] 

and 

/([4.5,5]) = [20.25,25] - [18,20] -b 3 = [3.25,10] 
and we have our resullQ. 

2.2. Description of algorithm. We first note that we need not search in all of 
CT > 1 to find zeroes of Cn{s). Spira [13] Thm. 1] proved that all zeroes of C,n{s) 
must have real part less than 1.85; this was sharpened in [5] Theorem 3.1] to 1.73. 
We therefore need only consider cr G (1,1.73]. We can improve this for some values 
of A, but, as we shall see in 112.31 this is more than sufficient for our purposes. 

Let us consider the case N = 28. Let p denote a prime and let Op = tlogp. 
Hence we have 

exp(-f 6 » 2 ) , exp(-f 6 i 3 ) exp(-f 26 i 2 ) , , exp(-f( 26»2 -b 6 » 7 )) 

C(cr + tt) = l + - — - -b - — - -b - — - -b . . . -b - — - 

^Note that if we had written f{x) = [x — l)(fr — 3) then /[4, 5] = [3,4] • [1,2] = [3,8] which is 
the “correct” result. This sensitivity is common in expressions involving intervals. 
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and we will now write ( 28 ( 0 ', O 2 , ■ ■ ■ ,^ 23 ) for C 28 (o' + it) under such a change of 
variables. 

It would appear that we need to examine the space a G (1,1.73], 9p G [0, 27r) for 
p < 23, for zeroes. In fact we can do considerably better. First, we observe that 
On, dig and dgs only appear once in the sum. Call the sum without those three 
terms ^ 28 ' ■ ■ ■, O 13 ). Then C 28 cannot have a zero if there is no 0 , 62 , ■■■, di 3 

such that 

|C28'(lT,d2,...,di3)| < 17 -" + 19 -" + 23 -". 

We can go further. The dn and dig terms only appear on their own or in 
conjunction with d 2 . We write a = 11““^, b = 22“'^, c = 13“'^ and d = 26““^. Then 
a little high school geometry (the cosine rule to be precise) tells us that 


exp(-idii) exp(-i(dii + d 2 )) 
11 '^ 22 -'^ 


< \/a^ + 52 + 2 ab cos 02 , 


and 


exp(-idi 3 ) exp(-i(di 3 + d 2 )) 


< \/+ 2 cd cos 92 - 


13 '^ 26 -'^ 

Call C 28 " (o’, 92 , Oz, dg, Oj) the result obtained by removing the n = 11,13, 22 and 
26 terms from (( 28 '■ With a, b, c and d as above, define 


f{a, 02 ) = l7 + 19 '^ + 23 + 6 ^ + 2ab cosd 2 + + 2cd cos 92 - 


Then ((28 cannot have a zero if there is no ct, d 2 ,..., dy such that 

|C28"(o’, 02,^3 j^*5j^ 7)I </(0’,d2). 


Our algorithm is as follows. Divide a, 02 , Oz, 9z and 9^ into small intervals that 
cover [1,1.73] and [0, 277 ]"^ respectively. We refer to any choice of five such intervals 
as a “box”. Push all possible boxes onto the stack. While the stack is not empty, 
pop off a box and compute an interval z containing ]C 28 "| for that box. Compute 
an interval containing /(ct, d 2 ). If the interval z — /(ct, d 2 ) is wholly positive, then 
that box did not contain any zeroes, so discard it. If the interval is wholly negative, 
then terminate with failureo If the interval straddles zero, then divide the box into 
16 smaller boxes by halving the intervals for the dp, and push these new boxes onto 
the stack. 


2.3. Details of the implementation. We implemented this algorithm in “C++” 
using our own double precision interval package written in assembler. This exploits 
an idea of Tambov [7] to make efficient use of the SSE instruction set of modern 
processors and uses CRMLIB m to implement the transcendental functions. 

We divided the interval for a into 16 sub-intervals [1 + (2" — 1) • 2“^®, 1 + (2"+^ — 
1) • 2“^®] for 0 < u < 15. Therefore, the first interval checked was tr = [1,1 + 2“^®] 
and the lasl[l was tr = — 2“^®, 2 — 2“^®]. Each of these intervals for a was handled 

by a single core of a compute node of the University of Bristol’s Bluecrystal Phase 
III cluster m Within a single core, the intervals for ^ 2 , ^ 3 , ^5 and 6 r were initially 


^We believe that this condition indicates the presence of infinitely many zeroes. We are grateful 
to the referee for suggesting a means by which one might seek to establish this, based on El,a 
and m- However, the weaker statement is sufficient for our purposes and we do not pursue this 
line of thought further. 

^Note that this covered a wider interval than was strictly necessary. 

single node of Phase III contains two 8-core Intel E5-2670 Sandybridge processors running 
at 2.6GHz. 
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divided into 16, 8, 4 and 2 sub-intervals respectively, for a total of 1,024 boxes. 
Since 02 contributed to more terms than the other variables, it made sense to start 
with a narrower search here: this seemed to work well in practice. 

Table [2] shows the data for N = 28 and a G [1,1 + 2“^®]. At each iteration, a 
box could result in 16 new boxes; at first this is what we see. We see that after the 
second iteration, the remaining search space decreases dramatically. 

Table 2. Number of boxes at each iteration for N = 28 and a G 
[I,l-k2-i6] 


Iteration 

Number of Boxes 

Coverage (%) 

1 

1,024 

100 

2 

16,256 

99.2 

3 

45,920 

17.5 

4 

118,560 

2.83 

5 

170,048 

0.25 

6 

195,920 

0.018 

7 

212,960 

0.0012 

8 

82,016 

0.000030 


We ran this algorithm for those N in (ED and in every case confirmed that Cn{s) 
has no zeroes for ct > 1. Checking each N took much less than a minute elapsed 
time using 16 cores, with = 21 taking the longest at 30 seconds. 
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